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Abstract 

We speculate about the possible interpretations of the recently observed excess in the h —» rp decay. We 
derive a robust lower bound on the Higgs boson coupling strength to a tau and a muon, even in presence of 
the most general new physics affecting other Higgs properties. Then we reevaluate complementary indirect 
constraints coming from low energy observables as well as from theoretical considerations. In particular, 
the tentative signal should lead to r —> py at rates which could be observed at Belle II. In turn we show 
that, barring fine-tuned cancellations, the effect can only be accommodated within models with an extended 
scalar sector. These general conclusions are demonstrated using a number of explicit new physics models. 
Finally we show how, given the h —» rp signal, the current and future searches for p —y ey and p —y e 
nuclear conversions unambiguously constrain the allowed rates for h —> re. 
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I. INTRODUCTION 


The discovery of the Higgs boson at the LHC [1, 2] imbues the standard model (SM) of par¬ 
ticle physics with completeness and self-consistency. Nonetheless, its failure to account for non¬ 
vanishing neutrino masses is one of the main motivations for considering physics beyond the 
SM. Incidentally, the accidental SM symmetries that prevent neutrinos from acquiring mass also 
completely suppress lepton flavor violating (LFV) processes. The observation of the former thus 
provides ample motivation for a rich experimental program to search for the latter. The CMS 
collaboration has recently reported a slight excess with a significance of 2.4 a in the search for 
LFV decay h —> r/i [3]. The best fit for the branching ratio of the Higgs boson to r/i final state 
(summed over r~/i + and t + /i - ), assuming SM Higgs production, is found to be 

B(/i ^ t//) = (0.84±°jj?) % . (1) 

This recent hint has expectedly received significant amount of attention in the literature [4-10]. 
It is thus an imperative to either confirm or reject validity of this tantalizing hint with more data 
by both ATLAS and CMS experiments. At the same time, it is instructive to revisit expectations 
for this observable within various new physics (NP) scenarios and in particular re-evaluate the 
feasibility of obtaining such a large signal in light of severe indirect constraints on LFV Higgs 
interactions coming from low energy probes. 

Without loss of generality, one can parameterize the mass terms and Higgs boson couplings of 
charged leptons after electroweak symmetry breaking (EWSB) as 

£y e = - mA/ 1 ' L f R - y i:j (t l L f R ) h + ... + h.c., (2) 

where the ellipsis denotes non-renormalizable interactions involving more than one Higgs boson 
and t l = e, //, r. In the SM, the Higgs couplings are diagonal and given by y tJ = {m r /v)8 lv 
where v = 246 GeV. On the other hand, non-zero y TjL and/or y /ir induce h —> t/i decays with a 
branching ratio of 

B(h r/i) = (W + \y^r\ 2 ) ■ (3) 

Assuming that the total Higgs boson decay width (TO is given by its SM value enlarged only 
by the contribution from h —> rn itself, i.e., = T/] M /[1 — B(h —> r/i)], where T^ M (m h = 

125 GeV) = 4.07 MeV [11], the measurement in Eq. (1) can be interpreted as a two-sided bound 
on the |i/r/J 2 +12/^r I “ combination of couplings. These limits read (see also the left-hand side panel 
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of Fig. 1) 


0.0019(0.0008) < yj \y Tf ,\ 2 + |y MT | 2 < 0.0032(0.0036) at 68% (95%) C.L.. (4) 

In general, the experimentally measured h -4 t/j event yield depends not only on the values of 
y Tll and y /iT , but also on other Higgs couplings contributing both to its total decay width 1% as 
well as its production cross-section (ah)- In particular, a given signal can be reproduced for larger 
(smaller) values of \y Tfl \ and y /lT by enhancing (suppressing) F/, and/or suppressing (enhancing) 
ah- Since both 1% and ay affect other currently measured Higgs observables, their individual 
effects can be disentangled by performing a global fit to all Higgs production and decay event 
yields at the LHC. The details of this procedure can be found in Appendix A, while the resulting 
best fit x 2 values as functions of \J\y T h\ 2 + \Vn rl 2 are shown in the right-hand side panel of Fig. 1 . 
In particular, ah is well determined by the measurements of both inclusive and separate exclusive 
Higgs production channels in several different decay modes. Similarly, T h is bounded from below 
by the observations (or at least indications [12]) of the dominant SM Higgs decay modes (h —* bb, 
h -x WW*, h —y r + r , etc.). Consequently the lower range of allowed y r/J and y //r values 
does not change much compared to Eq. (4) when considering the global fit. On the other hand, the 
fact that the total Higgs decay width is currently only weakly bounded from above [13, 14] allows 
significantly larger \y Tfi \ and y /IT couplings to reproduce the same observed signal in the general 
case. 1 Numerically, we find 

0.0017(0.0007) < yJ\y T h\ 2 + \y^\ 2 < 0.0036(0.0047) at 68% (95%) C.L.. (5) 

The rest of the paper is devoted to interpreting these ranges in terms of hypothetical new physics 
(NP) effects. In Sec. II we review the model independent considerations of flavor violating Higgs 
couplings using an effective field theory approach. Sec. Ill is devoted to a case study of the two 
Higgs doublet model (THDM) with a generic Yukawa structure (the so-called type-III THDM) that 
can account for the observed anomaly at tree level and be in agreement with present low energy 
flavor constraints. Scenario that relies on vector-like fermions to explain the data is presented in 
Sec. IV. Possibility to explain the anomaly via one-loop physics effects is discussed in Sec. V, 
where we investigate phenomenology of scalar leptoquarks to demonstrate difficulties with this 
1 It has been argued recently [15, 16], that the current total Higgs decay width measurements using off-peak data [13, 
14] introduce some model-specific assumptions into the fit. We have checked that removing the total Higgs width 
constraint from the fit, no upper bound on \f\y T h\ 2 + ly^rl 2 can be set. 
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FIG. 1. (Left-hand side panel) Higgs branching ratio to the r/i final state B(h —> t/i.) as a function of the 
effective /ir/r coupling. (Right-hand side panel) The x 2 fit to the latest available LHC Higgs data including 
the CMS h —> Tfi signal: varying only the effective hr/i coupling (shown in red, solid curve), and varying 
all Higgs boson couplings (shown in blue, dashed curve). 

type of approach. In particular, Sec. VD gives a very special solution that can accommodate 
the anomaly through an interplay between two different one-loop effects, and satisfy constraints 
imposed by low energy flavor physics experiments at the price of fine-tuning. We briefly conclude 
in Sec. VI. 


II. MODEL INDEPENDENT CONSIDERATIONS 

In the following we assume the SM contains all the relevant degrees of freedom at energies 
0(few 100) GeV. 2 In particular, any additional heavy field can be integrated out, so that we can 
apply effective field theory (EFT) methods. If the electroweak symmetry is realized linearly, 
then current data already strongly suggest that h for the most part should be a neutral CP even 
component of a linear combination of Y = 1/2 hypercharge weak doublets H a = (//,/, v a + x a h + 
.. ,) T , where the ellipsis denotes possible projections to other heavier scalar mass eigenstates and 
V/ ( \ x a\ 2 ~ 1/2. In general, EWSB can receive contributions from several sources (for example 
several Higgs multiplet condensates). However, electroweak precision tests severely constrain 
contributions from higher representations, so that also v 2 a < v 2 /2. 

2 In the minimal extension of the SM with Dirac neutrinos accounting for the observed neutrino masses, a non¬ 
vanishing h Tfj, rate is induced at the one-loop level. However, this contribution is negligibly small since the 
relevant decay amplitude is suppressed by two powers of neutrino Yukawa couplings in addition to the tau Yukawa. 
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In the decoupling limit we are considering, the contributions from the heavier Higgses are 
power suppressed and can be completely described in terms of higher dimensional operators. The 
leptonic Yukawa sector of such a theory, including the leading (dimension six) non-renormalizable 
operators suppressed by the EFT cut-off scale A, can be written as 

C Yt = - A'}/.,//„/•;, - A'f'+ h.c., (6) 

where L, and E r correspond to the leptonic weak doublets ( Li = (t L . z/j 7 ) and singlets (E, = t R ), 
respectively, and we explicitly show the flavor indices i,j = 1,2,3. Any additional dimension 
six operators coupling the Higgses to the leptons can be shown to be either redundant or not to 
contribute to the effective Lagrangian (2) after EWSB [17]. We can thus identify 

Vij @ij T ^ ij i (7) 


where 


e = V L 


\ a Va I ^ - 1 ) + I ^ + H + ±1-1 


A 2 


aVpVry 


Xr 


Xp 


Xr, 


Vp 


Vi, 


with v a = v a /v and the unitary matrices Vl,r diagonalize the leptonic mass term as 

™=Vl (^ a v a + A ,a ^v a vpv^j V R . 


( 8 ) 


(9) 


We can note immediately the two possible sources of non-vanishing y Ttl and/or y jlT . In a theory 
with multiple Higgs doublets, the first term in the square brackets in Eq. (8) can be non-vanishing 
if for at least two terms in the a-sum (with non-commuting A Q ) x a ^ v a . Thus (lepton) flavor 
violating Higgs boson interactions are possible even without additional NP contributions above 
the EFT cut-off. The simplest example of such a theory, the two Higgs doublet model (THDM) 
of type-III has recently been analyzed in this context in Refs. [6, 10] and we consider it in more 
detail in Sec. III. On the other hand, in a single Higgs theory only v\ = v/a/2 and x\ = l/\/2 
are non-zero so that the first term in the square brackets in Eq. (8) vanishes identically. Thus, 
the observation of /i — )■ r/i in this context necessarily implies the presence of the second term in 
Eq. (6). In particular, the CMS measurement can be interpreted directly in terms of the effective 
NP scale of 

0.84% ^ ' ' ll/4 


A ~ 4 TeV 


\B(h —> T/i) 


114 A 


+ \V L X 


/lllT/t|2 

V R\flT 


( 10 ) 


We observe that percent level B(h —> r/i) can be accommodated within the EFT approach (con¬ 
taining a single light Higgs doublet) with a sizable mass gap. At the same time, new degrees 
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FIG. 2. Comparison of the CMS implied range of hry couplings (green band) and the theoretical constraint 
imposed by the hierarchy of tau and muon masses (orange top-right region). See text for details. 

of freedom generating the observed h —> t/i excess should be within the kinematic reach of the 
second LHC run. 

Since A and A' contribute to both e as well as m the observed hierarchical structure of the 
charged lepton masses can be used to define the natural ranges of eij. In particular, without delicate 
cancellations between the various contributions in Eq. (9), the observed hierarchy between the 
muon and tau lepton masses implies [18, 19] 



(ID 


This constraint is compared to the CMS preferred range for \y Tfl \ and y jlT from Eq. (4) in Fig. 2. 
We observe that the two indications are not in sharp conflict and thus the observation of B(h — > t/i) 
at the percent level can in principle be explained by natural NP although a mild hierarchy between 
y Tfl and y /ir is preferred in this case. 

Another class of model-independent constraints comes from the fact that the dimension six 
operators contained in Cy f in Eq. (6) will in general mix with other operators contained in the 
effective SM. In particular, restricting the discussion to a single Higgs setup with H = Hi, 
X' will mix under charged lepton Yukawa renormalization into (II 1 Il f' affecting, for example, 
Higgs boson pair production at the LHC. However, as shown in Refs. [20, 21], such mixing is 
suppressed by three powers of the charged lepton Yukawa matrix and is thus (baring fine-tuned 
cancellations) completely negligible in practice. More phenomenologically relevant are the UV 
finite one- and two-loop contributions to the operators LH(a ■ B)E and Lr a H(c r • W a )E, where 
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<j iiv = i\ y M , y^j/2, /j /)jy and W“„ are the hypercharge and weak isospin field strengths, respectively, 
and r a are the Pauli matrices. These operators can mediate radiative flavor violating lepton decays 
as well as contribute to leptonic anomalous electric and magnetic dipole moments, all of which are 
already tightly constrained by experiment. The most stringent constraint comes from the r —>■ /xy 
decay [17], mediated by the effective Lagrangian 


£eE. = C L Q Ll + CrQr, y + ll.C. , 


( 12 ) 


where Ql. Rl = (e/87T 2 ) r m T (fia a ^P Lt RT)F a /3, P L ^ R = (1 =F 7s)/2 and F a p is the electromag¬ 
netic field strength tensor. The Wilson coefficients cl,r receive comparable one- and two-loop 
contributions. In the experimentally justified approximation y /7// <7 y TT (both assumed real) and 
m /t <C m T -C m h , they are given by [17, 22, 23] 


(1—loop) 


^(2—loop) 


- Z2 Vr^y 


m 


1 1, 

-X + ylog — 
3 4 m 


ml 

2 


^(1 loop) 
R 


mi 


mi 


-0+7 log =2 


(125 GeV) 2 ^ 


y* (0.11 - 0.082^) , c 


(2—loop) 
R 


(125 GeV) 2 


m% f 

(13) 

Vnr{ 0-11 - 0.082 y u ) , 

(14) 


where y tt is the top quark Yukawa with the SM value of y tt = m t /v = 0.67 for a MS top mass of 
rh t — 164 GeV. The resulting EFT correlation between B(h — y Tfi) in Eq. (3), and B(t —y /n) 
given by 

K( \ T r a EM m r /< \2 , \ |2\ 

B[r ->■ fin/) = 64?r4 (\c L \ + \c R \ ) , ( 15 ) 

is shown in Fig. 3 (diagonal dashed orange line), assuming SM values of all Higgs boson couplings 
except y Tfl and y //r . In the same plot, the CMS preferred range of B(h -f t/i) in Eq. (1) is displayed 
by the horizontal blue band, while the current (£>(r —> /ry) < 4.4 x 10 -8 @ 90% C.L.) [24] and 
projected future ( B(t —r /xy) < 3 x 10 -9 @ 90% C.L.) [25] indirect constraints are shaded in 
light and dark gray vertical bands, respectively. We observe that within the EFT approach, the 
CMS signal is well compatible with the non-observation of r —» /xy at the B factories and will 
marginally remain so even at Belle II. This is in contrast with the situation in most explicit NP 
models generating non-zero y T/ , or y /lT , as we demonstrate in Secs. III-V . 

Before going through explicit examples, let us consider the appearance of (lepton) flavor vi¬ 
olating Higgs interactions from a symmetry point of view. In absence of Higgs Yukawa inter¬ 
actions, the SM possesses a large global flavor symmetry Q F = SU(3) 5 (where we have omit¬ 
ted the f/(l) factors). In particular, the leptons transform as L (3,1), E rs_/ (1,3) under 
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FIG. 3. Correlation between B{h —> t/i) and B(t —> /jry) in various NP scenarios. The present experimental 
result for B{h —» t/i ) is shown in horizontal blue band [3]. Current and future projections for B(t —y fry) 
experimental sensitivity are represented with vertical light [24] and dark [25] gray bands, respectively. 
Superimposed are the predictions within the EFT approach (diagonal dashed orange line), in the type-III 
THDM (green and black bands), in models with vector-like leptons (diagonal dotted purple line) and in 
models with scalar leptoquarks (diagonal red and orange shaded band). See text for details. 


Qp = SU( 3)l x SU(3) e G Gf- In the SM (without neutrino masses), the charged lepton Yukawa 
matrix A ~ (3, 3) is the only source of Qp breaking. Consequently all lepton interactions are 
flavor conserving in the charged lepton mass basis. Conversely, as also demonstrated explicitly 
in Eq. (8), the generation of lepton flavor violating Higgs interactions requires at least two non- 
aligned sources of lepton flavor symmetry breaking. At the tree level, there are only two possi¬ 
bilities: (1) one can enlarge the SM scalar sector, such that more than one Higgs doublet couples 
to the leptons (corresponding to the first term in Eq. (8)); (2) one can extend the leptonic sector 
by vector-like fermions, whose Dirac masses and mixing terms with SM chiral fields can pro¬ 
vide additional sources of Qp breaking. This leads to the appearance of the A' contributions after 
integrating out the new heavy fermionic states. Both possibilities are explored in the following 
sections. Example of an enlarged Higgs sector is given in Sec. Ill whereas the vector-like fermion 
case is discussed in Sec. IV. 


8 










Of course many more possibilities exist at the loop level (e.g. [26, 27]). In practice however, all 
such models predicting sizable LFV Higgs interactions necessarily suffer from a severe fine-tuning 
problem, if they are to simultaneously avoid the stringent r —> /xy constraint. We demonstrate this 
on an explicit example in Sec. VD. A general heuristic argument however goes as follows. The 
operator LHE(WH) has exactly the same transformation properties under both Qg and chiral 
transformations of leptons as operators LH(a ■ B)E and Lr a H(a ■ W a )E. Thus, the simultaneous 
presence of the former and absence of the latter cannot be protected by a symmetry valid in the low 
energy theory. If LHE(H^H) is generated at the loop level, it will necessarily involve charged 
states propagating in the loop. These charged states will couple to photons generating in turn c LiR 
in Eq. (12) at the same loop level as y r/ , and y /iT , respectively. The two sets of matching amplitudes 
will also exhibit identical dependence on the fundamental flavor parameters of the model. The 
resulting expectation cl,r ~ %"!J T ,i,,i T /vrn T is clearly in conflict with existing experimental results, 
since a percent level B{h —>• r/x) would imply even an order of magnitude bigger B(r —> /xy)! 
The accidental cancelation required in the matching procedure for the radiative decay in order to 
accommodate both results is thus expected to be approximately one part in 10 3 at the amplitude 
level. As we demonstrate in Secs. Ill, IV and V, these expectations are confirmed in explicit 
models. 


A. 


h —> r/x vs. h 


re 


A positive experimental indication for the h r/x decay when combined with existing strin¬ 
gent experimental limits on/x-e LFV processes can be used to constrain LFV r — e processes. In 
particular, the product of the h —> re and h —> r/x branching fractions is constrained from above 
by the rates of /x —» ey, and /x — e conversion on nuclei. This is due to the fact that tree level Higgs 
decays to r/x (re) depend on (y er ,re) while the same sets of couplings contribute at the loop 
level to /x —> ey and /x — e conversion via diagrams with a virtual r. 

The contributions to the /x —>• ey process stemming from a virtual r are m T enhanced with 
respect to diagrams with intermediate /x or e states [17] leading to 


r >' 

C L ~ 


m T —3 + Ax r — x^. — 2 log x 7 


y^Vre > 


mz 


x 7 = 


2 ’ 


( 16 ) 


7Tin 8m|(l — x T ) 3 

where we have neglected the effects of the light lepton masses while m T dependence is retained in 
its exact form. The coefficient C R is obtained from Eq. (16) by replacing y KI — y y* t . The /x —> ey 
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branching fraction is thus sensitive to a distinct combination of the LFV Yukawas: 


Bill -+ e 7 ) ~ B^ {ly^Vrel 2 + \y^y e r\ 2 ) , B= 185 . (17) 


On the other hand, y — e conversion on nuclei is most sensitive to vector current effective operators 
(e^vPL’Rfi) (qYq) which probe different combinations of LFV Yukawas. The largest contribution 
to the vector coupling of the proton follows from Eq. (A17) of Ref. [17] and reads, in the limit of 
massless light leptons, 

„( p ) a -16 + A5x t - 36xl + 7x 3 + 6(— 2 + 3x T ) logx r * /10 . 

& - X,nl( 1 - * T )* -' US) 

The above effective coupling is evaluated in the q 2 —» 0 limit, where q is the momentum exchanged 

(v) 

in the conversion process. As before the right-handed current coefficient (j-fiv ' s obtained from the 
above result by replacing y %3 with y* % . Note that due to strong experimental limits on B{y —* e 7 ) 
we can at present safely neglect the dipole operator Ql,r contributions to y — e conversion [17]. 
The relative conversion rate B{ye) Au = r(/r —» e) A u/r ca pture a u is then given by 


Biy -+ e) Au ~ 


(i++ + isSi 2 


rcapture Au 

= K (If/er^rl 2 + bre?/r M P) , 


B£ e = 4.67 X 10" 4 , 


(19) 


where the coefficient V <p) = 0.0974 m b J 2 and the muon capture rate r eaptureAu = 13.07 x 10 6 s _1 
have been taken from Ref. [28]. 

The complementary information on the LFV couplings extracted from the rates of y —> e 7 and 
y — e conversion is sufficient to completely correlate the Higgs LFV h -A ry and h -y re decays: 


Bih —» ry) x B(h —> re) = 


m h 


-i 2 


8tt r h _ 

= 7.95 x 10' 10 


Biy -»• e 7 ) Biy -+ e) Au 


B^n 

Bjy -A e 7 ) 

io- 13 




+ 3.15 x 10 


-4 


Biy —> e) Au 

lfF 13 


( 20 ) 


where in the second line we have approximated ~ r| M . The best experimental limit on B(qi —> 
e) Au < 7 x 10 -13 (at 90% C.L.) was achieved by the SINDRUM II Collaboration [29] while the 
best upper bound on Biy —> e 7 ) < 5.7 x 10~ 13 (at 90% C.L.) was recently determined by the 
MEG Collaboration [30]. Note that with the current experimental data the sum on the right-hand 
side of Eq. (20) is completely saturated by the y — e conversion contribution. Combining the two 
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bounds with the central value for the h —» r/i branching fraction in Eq. (1) leads to an upper 
bound B(h —» re) < 0.26. This is above the current indirect constraint coming from searches 
for r —> e7 [24] which reads B(h —>• re) < 0.19 . In the near future the experimental sensitivity 


icr 2 

"S icr 4 

T 
-c 

10“ 6 

10“ 8 

10-8 10 -6 10 -4 10 -2 

FIG. 4. Allowed region in the B(h —> re)-B{h —> r/i) plane when experimental upper bounds on /i — > cy 
and /i — e conversion rates are taken into account. Pink region is permitted in the effective theory setting 
while the dashed line indicates how much the region will shrink if Mu2e and MEG II experiments see no 
signal events. Green region is allowed within type-III THDM model with rn,\ = 0.5 TeV and tan j3 = 10. 
Rulers indicate how much the region shrinks with increasing tan /3 or rn a , while dashed lines correspond 
to improved experimental upper bounds on /r —> e'y and fi — e as described in the text. Independently, an 
upper bound from B(t —y e'y ) applies on the B(h —y re) in the effective theory. 



to both /i —> e'y and especially /i — e conversion in nuclei is expected to improve significantly. 
In particular the DeeMe [31] experiment aims at a sensitivity of the order of 1CT 14 while the 
Mu2e [32] experiment could improve existing SINDRUM II bounds by four orders of magnitude. 
Given B(h —y r/i) at the percent level, such a measurement will indirectly probe B(h —* re) at 
the order of magnitude 10” 5 . 
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III. TYPE-III TWO HIGGS DOUBLET MODEL 


Several well motivated scenarios of physics beyond the SM (e.g. supersymmetric extensions 
or models addressing the strong CP problem via a PQ symmetry) require the addition of another 
SU(2)l doublet scalar to the minimal SM Higgs sector (c.f. [19] for a recent general review). 
The most extensively studied version of the THDM is the type-II model in which the absence of 
tree level flavor changing neutral currents (FCNCs) is imposed by the requirement that only one 
of the doublets couples to up-type quarks, while the other one couples to down-type quarks and 
charged leptons. FCNCs in the quark sector are relegated to the one-loop level via the interactions 
involving charged scalars, whose couplings to quarks follow the CKM structure of the SM. In 
this model LFV is absent as in the SM due to massless neutrinos and because only a single A" in 
Eq. (6) is non-vanishing. 

In order to accommodate LFV Higgs interactions one thus has to depart from this version and 
consider the THDM with a generic Yukawa structure (type-III THDM). In the following we restrict 
our discussion to the parameter space with a MSSM-like Higgs potential, for which the following 
relations hold [33] 

„ v u „ m 2 A + m \ 

tan/9 = —, tan 2 a = tan 2/3 —^^ , 

v d rn A - m z ( 21 ) 

2 2,2 2 2,2 2 
m H ±=m A + m w m H = m A + m z - m h , 

where 0 < /3 < tt/2 is the angle of the rotation matrix that diagonalizes the mass squared matrices 
of charged scalars and pseudoscalars and —tt/2 < a < 0 is an analogous angle for the neutral 
scalars. The vacuum expectation values of the doublets are denoted by v u and v d and vj\/ 2 = 
V v u + v/ ~ 174 GeV. In our input we set the mass of the Higgs boson to the observed value 
mh — 125 GeV and choose (3 and rriA to be free parameters of the scalar potential. Such a 
THDM can arise as a low energy effective theory of the MSSM in the limit in which all the heavy 
superpartners are integrated out. The LFV interactions of the Higgs boson then originate from 
diagrams involving heavy sleptons in the loop, c.f. [34, 35]. Note however that generically in the 
MSSM, LFV radiative lepton decay contributions will be generated at the same order as discussed 
in Sec. II, leading to tighter constraints, c.f. [36]. We thus do not attempt to interpret the type-III 
THDM model parameters in terms of the underlying MSSM parameters, but consider them as 
essentially uncorrelated, only subject to experimental bounds and theoretical constraints such as 
perturbativity conditions. 
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The relevant part of the Yukawa Lagrangian in the charged lepton mass eigenstate basis is [33] 

( 22 ) 


H k H + 

£ = + -^-H + ULj£R t i + h.c., 


where the Yukawa couplings can be written as 

.k m U 


y% k = - 4 ( x % tan ^ - x u ) • 


(23) 


The terms e < j l parametrize the off-diagonal charged lepton Yukawa couplings and are free param¬ 
eters of the model. Finally, the coefficients x k for H k = (II. h, A) are given by [33] 

x * = (— sin a , — cos a , i cos (3) , 

(24) 

x k d = (— cos a, sin a, i sin (3 ) . 

In addition to tree level exchanges of neutral Higgses we have loop contributions to FCNCs due 
to charged Higgs, whose interactions with leptons are parameterized by the following Yukawa 
couplings: 


yf 


v'-^in jr/.; 

3 = 1 


PMNS / m ii 


Vd 


3i '-ji 


4 tan p 


(25) 


A. 


TH 


Using preceding relations, at the tree-level the type-III THDM contributes to the effective Higgs 
couplings as 

e l 

VyuT (th) — (sin a tan j3 + cos a), (26) 

and we can translate Eq. (4) into alcr two-sided bound on the e £ parameters 

0.0027 < | sin a tan /3 + cos a\ I4J 2 + |e^ T | 2 < 0.0045 . (27) 


For reference we also give the expression for the branching fraction 


B[h —>■ r/x) = ——- (sin a tan /3 + cos a) 2 ( \e l \ 2 + |4j 2 ) • (28) 

lo7rl h 

We show in Fig. 5 the points in the plane spanned by (144 + I444 2 ar| 4 <n.A for which B(h —> 
r/i) = 0.84% when all other Higgs decay rates are held SM-like. Interestingly, the decay exhibits 
small dependence on /3 for large values of tan B. Indeed, an expansion around the decoupling limit 
in small £ = m| / m 2 A gives 


sin a tan + cos a ~ 


2 tan j3 (tan 2 /? — 1) 
(tan 2 (3 + l) 3 / 2 


V 


2m| 

m A 


(29) 
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For light pseudo-scalar masses (m^), mixing with the heavy scalar Higgs could affect the cou¬ 
plings of the light Higgs to the gauge bosons. However, the overall modification factor of hWW 
and hZZ couplings relative to the SM values is given by sin(/3 — a) and starts to differ from 1 
only at order 0(( 2 ): 

Ztlitg — l ) 2 9 

sin(/3 - a) ~ 1 - [tl + l)* * ' (30) 

We have checked that even for m,A = 150 GeV the hWW coupling deviates less that 2 % from its 
SM value, for any value of tan j3. We conclude that the effect of modified light Higgs couplings to 
gauge bosons can be safely neglected at the current precision of Higgs properties measurements. 
The more important effects of possible modifications of the light Higgs (tau) Yukawa couplings 
are discussed in the next subsection. 



FIG. 5. The values of the couplings e^ (e^ T ) as a function of tua for which B(h t/i) = 0.84% in 
type-III THDM, when all other Higgs decay rates are held SM-like. The decay exhibits small dependence 
on (3 for large values of tan j3. See the text for details. 


B. Constraints from r —> [.r/ 

The coefficients cl,r in Eq. (12) are generated via one-loop diagrams with charged or neutral 
Higgses. In the limit m 7/ = 0 (for the sake of consistency also yff* = = 0), these coefficients 
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are [33] 


^(1 -loop) 

■'R 


y 1 

^,3 48m ^ U 


mz 


y^Vrr* - y^yrri 9 + 6 lo g ~2^) 


m 


*2 j 


96m^ + 


J*.2 „ J 

e T/A/3 S /3 6 rr 


C 


(1-loop) ^ 
L — 


E 


48m? 7 o 

fe=l,2,3 l 


y%*y% 


9 + 6 iog 


m; 


?7l 


*2 j 



5 


(31) 


The suppression of the above result by two small Yukawa couplings is relaxed at the two-loop 
level, where Barr-Zee type diagrams require only a single LFV Yukawa coupling along with y tt . 
Therefore the contributions of the neutral Higgses at the two-loop level can dominate over the 
one-loop result. Here we employ the Barr-Zee contributions due to neutral Higgses calculated 
for n —» e 7 in Ref. [37] and adapted to r — > /xy in Appendix Bl. The 2-loop contributions to 
B(t —>■ n'y) involve exactly the same combination of e e elements as B{h —* r/x) in Eq. (28). 
For the charged Higgs contributions to the best of our knowledge no complete calculation exists 
in the literature. A partial result with the charged Higgs contribution inducing the H® 77 * or 
H^rjZ* vertices has been presented in Ref. [35]. However, the diagrams with an effective II = W T y 
coupling via loop of t and b are still missing. In the following we use the known partial results 
to estimate the order of magnitude of the constraint. The missing contributions could potentially 
further strengthen the r —> /ry bound but are not expected to substantially weaken it. 

In order to explore the nontrivial correlation between the B(h —* r/i) and B(t — » /j. 7 ) we 
sample the parameter space of , e^ r , e e TT for a specific choice of tan /3 and rn a ■ The ranges 
allowed for e^ /M , T are required to fulfill the naturalness criterium of Eq. (11) 

|e^ T e^| < ^ T (sin a tan j3 + cosa) -2 , (32) 

and that each of the effective couplings is within the perturbative range 


I e lir, rfi I < ■'/ 47 tv / 2 I sin a tan (3 + cos a \ 1 . (33) 

The 1-loop amplitude of the r —>■ /xy process further depends on the diagonal y^ T Yukawa 
coupling which is experimentally constrained by the searches for h tt decays within ATLAS 
and CMS experiments. In fact both collaborations report strong evidence for the existence of 
this mode with the signal strengths (normalized to SM expectations) ft TT = 1.43 l |j^'-i and /x rr = 
0.78 ± 0.27, observed by ATLAS [38] and CMS [39], respectively. In the following we use a naive 
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average of the two experimental results: 


/i TT — 1.02+qIo ■ (34) 

In the decoupling limit (c.f. Eq. (30)) and with //[) SM-like, we can assume SM-like Higgs pro¬ 
duction cross sections and model the departure of fi TT from 1 by an appropriate shift of y ] l T away 
from — a / 2 m T /v (i.e. non-vanishing e e TT ). 

A scenario with SM-like y h TT coupling corresponding to fi TT = 1 for fixed tan (3 = 10 and 
rriA = 0.3 TeV is presented in Fig. 3 by a black narrow stripe. This scenario easily passes both 
experimental constraints. A perfect one-to-one correspondence between the two observables is 
spoiled by the asymmetry between the left- and right-handed r —>■ /ry amplitudes at one-loop, see 
Eq. (31). On the other hand, for masses rriA significantly larger than 500 GeV it is not possible to 
reconcile both predictions with the corresponding experimental values. 

Allowing fJ TT to vary within the experimental bounds imposed by the h —> rr decays, i.e. 
Eq. (34), one obtains predictions for h —> r/r and r —>■ fi'y as represented by the green band 
in Fig. 3. In particular, this additional freedom in r —> fiy breaks, to some extent, the strict 
correlation with the /wr/i rate. The remaining correlation signals the dominance of the 2-loop 
contributions in B{r —> /n) which are independent of e e TT . 

The interaction vertices that induce the LFV decay of Higgs boson to Tfi pairs also contribute to 
the LFV decay of tau lepton to three muons, r~ —>■ /r /i + fir . The current best experimental upper 
limit on its branching fraction has been reported by the Belle Collaboration £>(V —>• fi~fi + fi~) < 
2.1 x 10~ 8 at the 90% C.L. [40]. This process proceeds through the tree level exchange of the 
neutral scalars /z, H and the pseudoscalar A. Since all the relevant amplitudes are suppressed 
by flavor-diagonal muon Yukawa, this constraint is not competitive with r —* /ry in terms of 
sensitivity to e £ rfl and 


C. h —)• Tfi vs. h —> re 

Here we repeat the derivation of the bound on h re described in Sec. IIA in the context of 
the type-III THDM model. In this case, both // —>■ ey and f i — e conversion receive contributions 
not from a single Higgs but from three neutral Higgses while the charged Higgs contributions 
are relatively suppressed by light lepton masses. The effective operator coefficients presented in 
Sec. IIA now have to be summed over the three Higgses with proper mass and effective coupling 
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replacements. Ultimately, the observables depend on tp = tan (3, rriA, and the LFV parameters 
ef ■, where dependence on the first two parameters can be absorbed into an overall coefficient: 


B{n ey) ~ B^ n (tp,m A ) (|e Mr e re | 2 + \e eT y Tfl \ 2 ) , 
—> e) A u — B% e (tp,m A ) (| Cer^i-iT | “i“ |^re^r/x| ) ? 


B^ n (t 0 = 10, m A = 0.5 TeV) = 2.67, 
H(( e (^ = 10, = 0.5 TeV) = 0.037. 

(35) 


Tree level LFV Higgs decays h —>■ r/i and h —* re are still properly described in the effective 
framework of Eq. (26). In contrast with what has been found in the effective framework (c.f. 
Eq. (19)), the sensitivity to ef ? - is now similar in both low energy observables. Also, the coefficients 
now have a strong dependence on tan (3 and m A . Finally, Higgs decays are suppressed 
by a projection factor so that 


B(h —* r/i) x B(h —* re) 


raft 
167rr ? 


(sin a tan /3 + cos a) 


/ B{fi ey) H(/i ->• e) Au A 
VH^ e7 (^,ra A ) B^ e (t 0 ,m A )) 

(36) 


These features suppress the h —>• re process very stringently in the type-III THDM model: 


B(h —> re) < 6 x 10 6 ; > 0.3 TeV , tan (3 > 2 . (37) 

The above bound is monotonically decreasing with tan (3 and r«/i as can be seen on Fig. 4. Note 
that at the projected sensitivity of the Mu2e experiment, the right hand side of Eq. (36) would 
start to be dominated by the /i —> ey constraint, so an improvement expected from the planned 
MEG upgrade (MEGII) [41] will be essential in this context, as can be seen by comparing the two 
bottom-most dashed contours in the plot. 


IV. VECTOR-LIKE LEPTONS 

In this section we discuss the possibility of generating sizable LFV Higgs couplings by mixing 
the SM chiral leptons with new vector-like leptons (VLLs). Such states can appear in the low 
energy spectra of grand unified theories [42-46] and are predicted to exist in composite Higgs 
scenarios with partial compositeness [47-50]. In this setting the chiral leptons obtain additional 
couplings to the Higgs by mixing with the heavy vector-like leptons. 

We first consider the addition of VLLs in a single vectorial representation of the SM gauge 
group SU(3) C x SU(2) l x U( 1) y , be. (1, 2)y 2 © (1,2) i /2 or (1, l)i © (1, l)_i. In this case one 
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can prove [51] that LFV Higgs couplings are directly related to LFV Z-boson couplings 


r z - 

Z -'LFV 


fj 


2 cw 


(W^ - Vkr^i) z „, 


( 38 ) 


where g ~ 0.65 and cw = cos 9 W ~ 0.88 are the weak isospin gauge coupling and the cosine of 
the weak mixing angle, respectively. In particular one obtains a relation = ( rrij/v)(Xij — Y l;J ). 
Severe constraints on V T/i/iT , Y TI , flT < 1CT 3 from searches for r —y 3/i [52] then preclude any 
sizable contributions to B(h —> t/Y) ? De-correlating LFV Higgs and Z boson couplings requires 
the introduction of VLLs in both SM gauge representations ('F B , V F L ) mixing with chiral leptons 


Cvll = X^ e H{1 - 75 )^ l + X^ E H{1 + 7s )T L 

+M* (A e E^ E + XiL^ L + C L V L ty L + C r ^ e ^ e ) + h.c.. 


(39) 


Integrating out the vector-l ik e states one finds the LFV Higgs and Z boson couplings can now be 
completely de-correlated in the limit where the direct couplings of chiral fermions to the Higgs 
vanish, and the SM leptons obtain their masses exclusively through mixing with VLLs in Eq. (39). 
In particular the lepton flavor off-diagonal Higgs Yukawas in Eq. (7) are then of the form [54] 

8v 2 


e = 


XiC^XvC^XvC^XvC-^X 


■Y-l 


-1 • 


1-1 • 


M 2 


(40) 


All the LFV phenomena are driven by this contribution and one can derive a one-to-one correlation 
between tree-level Higgs decay h —)■ tji and the 1-loop radiative decay r —>■ /i 7 [54] 

B(h —> r/i) 47 t B{h — > t + t~) sm 


2 x 10" 


(41) 


B[t —>■ /i 7 ) 3a B[t —>■ /iFi/)sm 

The 1-loop suppression factor of r —)• /i 7 is not small enough to evade the experimental upper 
bound and at the same time accommodate the h —> r/i at the percent level, as illustrated by a 
pink-dotted line in Fig. 3. 


V. SCALAR LEPTOQUARKS 

A scalar leptoquark state (LQ) can induce h —> r/i decay via quark-LQ penguin diagrams. 
This type of diagrams requires a helicity flip on one of the fermion lines. One therefore expects 
suppressed amplitude when all the fermions are substantially lighter than v. It turns out that size 


3 A similar conclusion can be drawn in the case of weak-isospin triplet Majorana fermions, appearing in Type-III 
see-saw models of neutrino masses [53], 
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FIG. 6. Feynman diagrams for LQ contribution to h — > t/i. 


of LQ Yukawa couplings must then be well beyond the perturbative regime to result in B(h —> t/j ,) 
at the percent level. In this respect, only LQ states that couple to charged leptons and top quark 
are suitable candidates that will be studied in what follows. Inspection of helicity structure of the 
diagrams reveals that both left and right chiralities of leptons and top quark have to couple to the 
LQ state. 

The physical Higgs can couple to the scalar A or to the top quark as shown in the first two 
diagrams in Fig. 6 . While the strength of the coupling relevant for the latter process is fixed by 
the top Yukawa, the former process depends on an unknown (iLQLQ coupling, Xv, that originates 
from the marginal “Higgs portal” operator, 

C 3 -XH ] HA ] A . (42) 

Here A is the scalar LQ and H is the SM Higgs doublet. 


A. The Ai = (3,1, —1/3) case 

The Yukawa couplings of Ai are given by the following Lagrangian 

£ Al = yijQ i,a A±e ab L c i' b + y^U i AiE Cj + h.c., (43) 

where Q l = (u' L , d l L ) T and U l = u' R are the quark weak doublets and up-type singlets, respectively. 
We explicitly show flavor indices i, j = 1, 2, 3, and SU{2 ) indices a, b = 1, 2, with e ±2 = 1. Also, 
here y/j and y T X are elements of arbitrary complex 3x3 Yukawa coupling matrices. After expanding 
the SU (2) indices, we obtain 

= y-^VA, - {V^ KM y L V VMm ) i3 d' L v c L >Ai + y'^A'/A: + h.c., (44) 

where Vckm and Vpmns represent Cabibbo-Kobayashi-Maskawa and Pontecorvo-Maki-Nakagawa- 
Sakata mixing matrices, respectively. All fields in Eq. (44) are specified in the mass eigenstate 
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basis. The Wilson coefficients of the effective Lagrangian (3) are obtained after summing 

the diagrams in Fig. 6 

2 /t M (jit) = ~rn Al ) y*yg (y£y £) , (45) 

where N c = 3 is the number of colors. The relevant loop function further depends on the portal 
coupling A, 


gi{X,m A ) = (m 2 A +m 2 )C o (0,0,m 2 h ,m 2 ,m A ,m 2 ) + B 0 (m 2 h , m 2 , m 2 ) - B 0 (0,m 2 , m] 


(46) 


+ An 2 C 0 (0, 0, ml, m 2 A , m 2 , m A ), 

given in terms of Passarino-Veltman functions B 0 and 6 ' (l 4 . The h —» r/i decay width in the A, 
LQ scenario is then given by 

Qm^m 2 


T(h —y Tji) = 


\ytvtr\ 2 


+ \yirvl\ 2 ) |0i(A,m Al )| : 


(47) 


2 13 7T 5 U 2 

The state Ai also contributes to the r —^ 7 through leptoquark Yukawas present in h —> r//. 
The Wilson coefficients of effective Lagrangian introduced in Eq. (12) are 


C L(R) - 


Nc mt , . R* L* / R L \ 

hl ^ y ^ VtT ( y tT y tn), 


h(x ) = 


7 — 8 x + x 2 + (4 + 2x ) In x 
(1 — x ) 3 


(48) 


with x t = m 2 /m 2 Ai , in agreement with [56]. The r —* /17 branching fraction in presence of Ai 
depends on identical Yukawa combination as the h —>■ tji decay rate: 


(3(r /X7) = 


am; 


mi 


2 12 7r 4 T r 


^i(^) 2 (|s/&2/£| 2 + \ytyZ\ 2 ) 


(49) 


Additional constraints on the model parameters can also come from the measurements of (g — 
2)^ and the electroweak precision observables, i.e. B(Z —y bb). The Ai contribution to the former 
is dominated by the coupling combination which can always be suppressed. For the Z —>• bb 
the presence of Ai will modify the effective coupling Sg h L , defined as: 


£-Zbb ~ 


COS 9 W 


Z“H H + Sai)PL + (A + ^r)Pr] b 


(50) 


Here g is the SU (2 ) L gauge coupling, 6 W is the Weinberg angle. g b L R are the SM effective cou¬ 
plings including the electroweak radiative corrections, while 5g b L R parameterize contributions that 
originate from beyond the SM. The LQ state Ai contributes to 5g b L through loops with neutrinos, 


Sg b L = 


m 2 (is/fti 2 + k\ 2 ) 


288ir 2 


3 log x + 1 + 3 * 7 t sin 2 6w 


x 


3x 


( 51 ) 


4 We disagree with the result of Ref. [55] where the authors employed a subtraction at kinematical point mu = 0 to 
render the h —>• rg amplitude finite. 
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where x = m 2 Ai /m 2 z 1. The constraint on LQ Yukawas follows from combining Eq. (51) with 
the results of the recent fit to the Z-pole observables [51]: 

R e[5g b L ] = (2 ± 2 ) x 1 (T 3 , Re[<5^] = (15 ± 6 ) x 10~ 3 . (52) 

The Z- pole measurements are sensitive to interference terms between g b L R and 5g b L R , where g b L R 
are approximately real numbers, and therefore only real parts of 8g b L R give linear shift in g b L R . 
According to Ref. [51] the SM point corresponds to ^sm = 7.5, where 0 is defined as the fit 
with free 5g b L R . The presence of LQ will push 5g b L in the negative direction. At lcr we insist on 
X 2 < Xsm + 1 that corresponds to a lower limit A g b L > —1.2 x 10~ 4 . An approximate constraint 
on the Yukawas is finally obtained, 

^J\Vtr\ 2 + K\ 2 £ 1.30 ^ + 0.34 , (valid for 0.5 TeV < m Al < 2.0 TeV). (53) 

Clearly, this constraint does not preclude the Yukawa couplings combination entering B(h —» r/x) 
from being large. 

In the left-hand side panel of Fig. 7 we show the dependence of B{h —» r//) on the leptoquark 
Yukawa couplings assuming the total Higgs decay width to be SM-like. Here, we turn off the 
Higgs portal coupling. We set the mass of LQ to 650 GeV, consistent with current direct search 
limits at the LHC [57, 58]. In the same panel of Fig. 7, we show the CMS measurement reported 
in [3]. The B(t —> //y) dependence on the leptoquark Yukawa couplings for Ai case is shown in 
the right-hand side panel of Fig. 7 (dashed line), again for rri A] = 650 GeV. 

B. Impact of the Higgs portal and its constraints 

Next we study the impact of a non-zero Higgs portal coupling in the scenario with A^ As an 
example, for m Al = 650 GeV the loop function dependence on the portal coupling is 

gi(X, m Al = 650 GeV) = -(0.26 + 0.12A). (54) 

Thus, a positive large A could in principle relax the leptoquark Yukawa couplings and yield sizable 
h —» r/x rates without violating the r —>■ gq constraint. However, the Higgs portal coupling 
also induces corrections to the h —> 77 decay and to gluon-gluon fusion (ggF) induced Higgs 
production with the leptoquark running in the triangular loop. The modified ggF production cross 
section normalized to its SM value is 

= |c 9 | 2 , where c g = 1 + 0.24-^-V A iC'(r A ). (55) 

a ggF m A 
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FIG. 7. (Left-hand side panel) Higgs branching ratio to r/r final state in the presence of scalar leptoquarks 
Ai = (3,1, —1/3) and Ao = (3, 2, 7/6) with Higgs portal coupling A set to zero. (Right-hand side panel) 
r —> q 7 constraints. The Ai (A 2 ) case is rendered in dashed (solid) line. The following transformation 
needs to be applied when going from the Ai to the A 2 case: y/j —y y^. 


Here, N A , is the number of A* components in the weak multiplet A. C(r A ) is the index of color 
representation r A of A and for the triplet ((7(3) = 1/2). We consider heavy enough colored 
scalars such that the loop function is in the decoupling limit. Similarly, the modified h —> 77 
decay width, normalized to its SM value, is given by 


f h—t 77 
■pSAf 
h —^77 


= C~ 


Xv 2 

where cL = 1 — 0.025— ^-d(r A )'S~'Q 

mi 


(56) 


The sum in Eq. (56) runs over all weak components of the SU(2) L multiplet. d(r A ) and Q A i are 
the dimension of the color representation of A and the electric charges of weak A* components, 
respectively. We fit the latest LHC Higgs data (including the CMS signal of h —y Tfi) taking A 
and leptoquark Yukawa couplings as free parameters. The result for the Ai LQ model is shown in 
Fig. 8 . The preferred regions at 1 a and 2 a are rendered in solid pink and dashed pink, respectively. 
If besides the A and leptoquark Yukawa couplings, we also allow other Higgs couplings to vary 
(see Appendix A) we instead get the dark and light grey regions at 1 a and 2 a levels, respectively. 

The correlation between the h —> T[i and r —> /ry branching ratios in presence of the 
(3,1,—1/3) leptoquark state for |A| < 1 and m Al > 600 GeV is depicted in Fig. 3 with a 
pink stripe. 
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FIG. 8. Fit to the latest available LHC Higgs data including the CMS h t/j, measurement. 


C. The A 2 = (3, 2 . 7/6) case 

The Yukawa couplings of the A 2 leptoquark to SM fermions are 

£a 2 = yl-E'Arcr - yWA a 2 e ab L^ b + h.c., 


(57) 


71J Z '"V UlJ 

where we explicitly show the flavor indices i, j = 1, 2, 3, and SU(2) indices a, b = 1, 2. y L and 
y R in Eq. (57) are arbitrary complex 3x3 Yukawa matrices. In the mass eigenstate basis we have 

£a 3 = <444Af - + (FV* KM ) y 44Af- + (y R V'pMN S )«fi , R 4Af - + h.c., 

(58) 

where a superscript on A 2 denotes the electric charge of a given SU (2) doublet component. Note 
that only the A^ 3 state couples to charged leptons and the top quark and thus contributes to 
h — > t /i decay via virtual top quark. Due to the weak-doublet nature, this state has an additional 
Higgs portal operator of the form (7/7 A) (ZJII) that does not induce A) /,! coupling to the physical 
Higgs field. The effective hry couplings are 

N c m t 


Vt/i At) 


m A 2 )y^tUtr (y%y%?) , 


(59) 


167T 2 V 

where the loop function gi has been introduced in the previous section. The h — > r/r decay rate is 
then 

r {h ->• r/i) = ^^| |gi(A,m A2 )l 2 {\y%ytr? + lz/rt?/£l 2 ) • ( 6 °) 


23 












Allowed values for the Higgs portal coupling A can be inferred from a global fit to the Higgs data 
as has been done for the portal coupling of the (3,1, —1/3) state. (See Fig. 8.) We do not attempt 
to repeat the same procedure for the state (3, 2, 7/6) since we do not expect a drastic change in the 
allowed range of A. 

The Wilson coefficients for r —> 7x7 are again proportional to the couplings responsible for 
h —* t/j: 


C L(R) = 


N c m t 
24 m T m\ 2 


h 2 (x t ) (y^*yrt) , h 2 (x) = 


1 — 8x + 7x 2 + (4 — ICte) In x 


(1 -xf 


( 61 ) 


with xt = m 2 /m 2 ^ and agree with the formulas presented in [56]. Finally, the corresponding 
branching fraction is given by 


B(t ->• fj.'y) = 


ami 


mi 


2 12 vr 4 r r 


h 2 (x t ) 2 {\y? T y^t\ 2 + 


(62) 


As in the case of the Ai leptoquark, contributions to the muon anomalous magnetic moment 
(g — 2) /; are proportional to y^y^ t and can be always suppressed. On the other hand, contributions 
of A 2 enter the Z —» bb via positive contributions to 5g b L that can actually help in easing the 
tension with the experimental data [59]. Most importantly, these observables do not constrain 
large possible Ai effects in h —>■ t/x . 

In the left-hand (right-hand) side panel of Fig. 7 we show the B(h —> t/j) (B(t —> 7 x 7 )) de¬ 
pendence on the leptoquark Yukawa couplings for the A 2 LQ case with the Higgs portal coupling 
turned off and taking ma, = 650 GeV. Also in this leptoquark scenario the bound on B(t —>• 7x7) 
excludes sizable B(h —> r/x) due to the strict correlation between the two observables. See the 
orange stripe in Fig. 3, where the portal coupling is restricted to |A| < 1. 


D. Fine-tuning solution 

In this section we give an example of a phenomenologically viable model with fine-tuning. 
We add to the SM a scalar leptoquark A x = (3,1, —1/3) and a vector-like top quark partner 
T' l , T' r = (3,1, 2/3). In the following, the top partner will mix with the SM top quark. In full 
generality, the relevant Lagrangian is given by [51] 

-CD y t q' 3L Ht' R + y T q' 3L HT R + M t T' l T r + h.c. . (63) 

After EWSB, the mass eigenstate fields t and T are obtained after rotating t' and T' with left- and 
right-handed mixing angles, 0 L and 0 R . In the phenomenologically viable limit of a heavy top 
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FIG. 9. Fine-tuning r —> py to zero while keeping h -X t/i at best fit point within the LQ model with a 
vector-like top quark partner. See text for details. 


partner (and small mixing), the following relations hold [51] 

m t ~ y t vl\pl , rn T ~ M T , sin 0 L « , sin^R ~ -^sin6* L . (64) 

m T t/ t ttiT 

More precisely, the direct LHC searches for a singlet top partner set a lower bound on rriT, while 
the electroweak precision observable (p parameter) constraints the mixing angle. In what follows, 
we set rriT = 700 GeV and sin 6^ = 0.2 which are the borderline values compatible with the 
constraints [51]. The relevant LQ interactions are 

r 3 vlfl'l L ^L c >- b + + 4j f ' R ^E Ci + h.c.. (65) 

After rotating to the mass basis 


C D y^-(cos6» L t L + sin6> L T L )if • 7 Ai-(-2/^(cos6> fl t fi + sin6» fl T fl )t'^' , Ai+x^(cos6» fl T R -sin6> fl t fi )£^ J Ai+h.c.. 

( 66 ) 


Now, r —>• p 7 amplitude will receive contributions both from t and T at one loop level. Assuming 
non-zero y 33 , y 32 and x§ 2 , neglecting the right-handed mixing and linearising in sin 6l, we find the 
branching ratio 


Q!EM® 7 1J/331 ^ 

e(T -*• ^ = 2^r r ml 


2/!<* m t hi{;rn 2 t /m 2 Ai ) + x 32 sin6» L m T hi(m 2 T /m^)) 2 , (67) 


where /?i(x) is defined in Eq. (48). As a numerical example, we take m,A 1 = 650 GeV, rrir = 
700 GeV and find the condition for the complete cancellation in B(t —>■ //yj to be y§2 = 
—0.63 x 32 sin 6^. 
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On the other hand, the interaction of the physical Higgs boson with the top quark and the top 
partner is given by 

— C D — cos 2 9 l t L t R h + —- sin 2 0 L T L T R h + —- sin 2Q L t L T R h + — sin 2 0 L T L t R h + h.c.. 

v v 2v 2v 

(68) 

We compute the decay rate for h —> r/i with t, T and particles running in the loops. After 
properly including all the diagrams exemplified in Fig. 6 (four vertex and four leg diagrams), we 
find a finite result. For the numerical benchmark, we get 

Uth ~ (°- 26y ^ + °- 43 x 32 sin Ol)v33 , ( 69 ) 

where we have linearised in sin 0 L . Finally, requiring the complete cancelation in r —» /iy (that is, 
yf 2 = —0.63 xf 2 sin 0 L ), we find = 0.47 gives the best fit to h ^ r/i excess. In passing we 

note that if rriT > rriA, the T —> A£ decays can produce spectacular signatures at the LHC. 

VI. CONCLUSIONS 

Prompted by the recent experimental hint of h —> r/i events by the CMS Collaboration, we 
have carefully examined the implications of LFV Higgs decays at the percent level on possible 
extensions of the SM. In particular we have shown how a tentative B(h —> r/i) signal can be 
combined with other Higgs measurements to yield a robust lower bound on the effective LFV 
Higgs Yukawa couplings to taus and muons. Then we have reexamined the connection between 
LFV Higgs decays and LFV radiative decays of charged leptons, and demonstrated using EFT 
methods that the current CMS hint implies r —)■ /iy at rates, which could be observable at the Belle 
II experiment. In explicit models, the r —>• /iy constraint is generically much more severe. In fact, 
an eventual observation of h r/i at the LHC together with existing indirect constraints would 
already single out an extended SM scalar sector as a required ingredient in any natural explanation, 
the minimal example being THDM of type III. We have also examined purely fermionic SM 
extensions and models where h —>■ r// is generated at loop level, only to show that without the 
introduction of extra Higgs doublets, reconciling all existing indirect constraints with percent level 
B(h —* r/i), when at all possible, requires a high degree of fine-tuning. Finally, we have shown 
how a positive signal of —)■ r/i can be combined with experimental searches for f t er decays 
and /i — e conversions in nuclei to yield robust bounds on B(h —> re). In particular, considering 
only low energy Higgs EFT effects, the two LFV Higgs decay rates could still be comparable. On 


26 



the other hand, the THDM III cannot accommodate both respective branching ratios at the percent 
level. Currently planned improvements in experimental searches for (i — e LFV processes will be 
able to probe the product B(h —> rfijBUi —> re) at the 10 7 level in generic EFT and to order 
10“ 12 or better within the THDM III. 
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Appendix A: Higgs fit 

In the following we explain the method used to fit the LHC Higgs data. Ideally, we want the 
measurements to be reported in terms of the signal strengths normalized to the SM predictions 

i _ °~( fc ) ^ (Alt 

P(fc) ,SM nSM ’ 
a (k) 

where index i represents a Higgs decay mode, while k denotes a production channel. Such ob¬ 
servables could then easily be expressed in terms of new physics parameters. However, the actual 
experimental categories are never pure and contain events from different production mechanisms. 
Furthermore, in order to fully reconstruct the total likelihood function, it is necessary to know 
all the correlations among the different categories, which are available only to the experimental 
collaborations. 

It has been argued recently in Ref. [60] that the existing measurements for a given decay chan¬ 
nel should be presented in terms of two-dimensional likelihoods, in which gluon-gluon fusion 
(ggF) and associated production with a top pair (ttH) are combined as one signal (i^(ggF+ttH)), 
while vector boson fusion (VBF) and associated production with a gauge boson (VH) as another, 
( mvBF+VH ))• This decomposition is particularly useful since, on one hand, ttH is sub-dominant 
with respect to ggF and poorly explored in the present data set, while VBF and VH receive com¬ 
mon corrections in wide class of NP models obeying the custodial invariance. Furthermore, the 
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TABLE I. ATLAS Higgs data as used in our fitting procedures. The separation into production mechanisms 
for a given decay channel is used if provided. 


Decay channel Production mode Signal strength Correlation & Reference 

ATLAS 

h->bb 

VH 

0.2 ±0.65 


[61] 

h -+ ZZ* 

ggF+ttH 

VBF+VH 

1.8 ±0.65 

-0.2 ±3.7 

P = 

-0.72, [62] 

h -+ WW* 

ggF+ttH 

VBF+VH 

0.82 ±0.37 

1.74 ±0.80 

P = 

-0.15, [63] 

h -+ 77 

ggF+ttH 

VBF+VH 

1.61 ±0.41 

1.87 ±0.80 

P = 

-0.28, [62] 

h —> TT 

ggF+ttH 

VBF+VH 

1.5 ± 1.6 

1.7 ±0.84 

P = 

-0.55, [63] 

h -+ invisible 

VH 

0.13 ±0.31 


[64] 

h -+ Z 7 

inclusive 

2.0 ±4.6 


[65] 

h -+ pp 

inclusive 

1.6 ±4.2 


[ 66 ] 


correlations are automatically provided by the experimental collaborations for a given decay chan¬ 
nel. 

Following the recommendation, ATLAS and CMS have combined different search categories 
for a given decay mode to provide separation into production mechanisms. The results are usually 
presented in 2D plots in (p( gg F+ttH), P-(vbf+vh)) plane. In this case, we parametrize the likelihood 
with 


Xi = 


E 


/ . \ T / . \ 

/ l ^ 1 \ / 1 S' 1 1 

| ggF+ttH) ~ P(ggF+ttH) | y-1 P(ggF+ttH) “ P(ggF+ttH ) | 


V A 4 l (VBF+VH) P'l VBF+VH) ) 


\ A 4 \vBF+VH) r(VBF+VH) ) 

where the sum goes over the decay channels and the covariance matrices are given by 

/ I ' 2 

Vi = 


(A2) 


P a [ggF+ttH) a (VBF+VH) 
2 


(A3) 


(ggF+ttH) y 

p l °(ggF+ttH)®(VBF+VH) y u (VBF+VH) 

We obtain the best-fit values (/)), variances (a) and correlations (p) from the plots provided by the 
experimental collaborations. These are presented in Tab. I and Tab. II. 
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TABLE II. CMS Higgs data as used in our fitting procedures. The separation into production mechanisms 
for a given decay channel is used if provided. 


Decay channel 

Production mode Signal strength Correlation & Reference 

CMS 


VH 

1.0 ±0.5 


[12] 

h -+ bb 

VBF 

0.7 ± 1.4 


[67] 


ttH 

1.0 ±2.0 


[68] 

h -A WW* 

ggF+ttH 

VBF+VH 

0.76 ± 0.23 

0.74 ±0.62 

P = 

-0.21, [69] 

h -> ZZ* 

ggF+ttH 

VBF+VH 

0.90 ±0.45 

1.7 ± 2.3 

P = 

-0.69, [70] 

h 77 

ggF+ttH 

VBF+VH 

0.50 ±0.41 

1.64 ±0.88 

P = 

-0.50, [71] 


0 -jet 

0.34 ± 1.09 



h —^ TT 

1 -jet 

2-jet (VBF-tag) 

1.07 ±0.46 

0.94 ±0.41 


[39] 


VH-tag 

-0.33 ± 1.02 




0 -jet 

0.77 ±0.55 



B(h -> T/r) [%] 

1 -jet 

0.59 ±0.62 


[3] 


2 -jet 

1.1 ±0.80 



h -» invisible 

VBF+VH 

0.14 ±0.22 


[72] 

h -A Z7 

inclusive 

0.0 ±4.8 


[73] 

h —> 

inclusive 

2.9 ±2.8 


[74] 


If the separation into production modes is not provided, we use the signal strength measure¬ 
ments from existing search categories. These in general target certain production mechanism, 
which, however, does not imply 100% purity. Inclusive categories are dominated by ggF (90%), 
while VBF-tagged categories can have 20% to 50% contamination from ggF. VH- and ttH-tagged 
categories are assumed to be pure. In this case, we write 


<?A-*h 


a 


SM 

A—th 


c a gg F 
^99^ SM 
a ggF 


+ ivBF 


&VBF 

SM 


a 


VBF 


+ ivH ~FT7 + ittH 


OttH 


O 


SM 

VH 


cr; 


SM ' 
ttH 


(A4) 
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where represent contributions of the specified production mechanisms to the given category. We 
do not consider correlations here and add each measurement to total x 2 separately, 



(A5) 


We approximate the likelihood for the total Higgs decay width measurement [75] by 



(A6) 


In Tab. I and Tab. II we finally summarize all the available LHC Higgs data used in our analyses. 

The total x 2 function is given by the sum of all the contributions, namely, ;\ /2 = xj + x! + Xr- In 
order to confront the new physics model to data, we express all signal strengths (//) and total Higgs 
decay width in terms of model parameters and minimize the x 2 10 find the best fit point. The best 
fit regions are defined by the appropriate cumulative distribution functions, namely, 68.2% (1 a) 
best-fit region for one- (two-) parameter fit satisfies x 2 — Xmin < 1-0 (2-3), while 95.5% (2cr) 
best-fit region satisfies 1.0 (2.3) < x 2 ~ Xmin < 4-0 (6.2). The other parameters in the likelihood 
are treated as nuisance parameters. 

The limits quoted in Eq. (4) are obtained after allowing only y Tfl , y llT to be free parameter while 
fixing other Higgs couplings to their SM values. On the other hand, in the limits shown in Eq. (5) 
we further allow for arbitrary values of the usual SM tree level Higgs boson couplings; K t , k t , n b , 
n v and k c where k v = n w = k z as well as arbitrary new contributions to loop induced Higgs 
couplings; k 9 and re 7 . 

Appendix B: Barr-Zee contributions to r —» /jy in type-III THDM 

The effective coefficients receive contributions due to the Barr-Zee diagrams 


Hrr(w) 
C R 


Hry(W) 


C R 


#77 (*) 



(B4) 


(B2) 


(B3) 


(Bl) 
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Here C H o ww measure the coupling strength of neutral Higgses to WW relative to the SM coupling 
hWW, i.e., deviation from the vertex igm w g Their values are C HW w = cos(/3 — a), C hWW = 
sin( / 5 — a). The coefficients shown above correspond to topologies that contain H k 77 coupling, 
mediated either by t or W loop. The virtual photon can be replaced by a Z boson, however such 
contributions are subdominant. Loop integrals /, g, and h were calculated in Ref. [37], whereas 
their arguments are defined as z t k = rn^/m 2 H o, and zwk = rri^/m 2 H0 . 
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